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Abstract. Metal-rich debris disks were discovered around several single DAZ and DBZ white 
dwarfs, which may stem from disruption of smaller rocky planetesimals. In some cases, the material 
in addition forms a gaseous disk. For SDSS J122859. 93+104032. 9, the double peaked infrared Ca II 
triplet at XX 8498,8542, 8662 A exhibits a strong red/violet asymmetry. With the Tubingen non- 
LTE Accretion Disk code AcDc we calculated the spectrum and vertical structure of the disk, 
assuming the chemical mixture of the disk's material being similar to a chondrite-like asteroid, the 
most prominent type in our own solar system. By modifying the code to simulate different non axis- 
symmetrical disk geometries, the first preliminary results are in good agreement with the observed 
asymmetric line profile. 
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MOTIVATION 
Metal-Rich Dust Disks 

Starting in 1987 with the analysis of G 29-38 by Zuckerman and Becklin [14], about 
twenty metal-enriched white dwarfs with a significant infrared excess in their spectra 
have been discovered, without a cool companion being found. Due to their high surface 
gravity the sedimentation time of these objects is rather short [8]. To explain the high 
metallicity by ongoing accretion, a metal-rich dust cloud was suggested as a donator 
source. Spitzer observations confirmed the dust material [12], which is located in an 
equatorial plane [4] forming a metal-rich but hydrogen- and helium-poor [6] dust disk. 



Metal-Rich Gas Disks 



For the hottest of those white dwarfs, with Teff > 20000 K, Kilic et al. [7] pro- 
claimed that their dust material is sublimated to a gaseous disk. And indeed, five 
white dwarfs, found in the Sloan Digital Sky Survey (e.g. Fig. 1) by Gansicke et al. 
[3, 2, 1, priv. comm.] indicate gas disk emission features in their spectra in addi- 
tion to the IR spectral component of a dust disk. These gas features are mainly the 
Ca II A A 8498, 8542, 8662 A, although at least two of the five spectra also show a Fe ii 
XX 5018, 5169 A in emission. 
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FIGURE 1. Spectram of SDSS J1228+1040 from the Sloan Digital Sky Survey. 

Origin of the Disk Material 

Former planetary systems around the white dwarfs progenitor may survive the host 
stars' late evolutionary phases, as different theoretical approaches have shown (e.g. Jura 
[5]). After the actual formation of the white dwarf, a smaller rocky body, an asteroid 
or planetesimal, might have been disturbed in its orbital motion by another surviving, 
larger planet. When it reaches the inner part of the system close to the central object, it 
will get tidally disrupted, providing the material to form the dust and gaseous disk. 

AIM AND METHOD 

Asymmetric Line Profiles 

The disk's spectral line features are expected to be double peaked as the emitting ma- 
terial rotates with its Kepler velocity around the white dwarf. Surprisingly, these double- 
peak structures show a significant difference in the red/violet peaks line strengths, 
shown in Fig. 2 for the prominent Ca ll feature of SDSS J122859.93-I-104032.9 (hereafter 
SDSS J1228-I-1040). We assumed this asymmetry to directly originate from an actual 
asymmetric disk geometry. Hydrodynamically simulating the shape of a disk changing 
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FIGURE 2. Can AA 8498,8542,8662 A in the spectrum of SDSS J1228+1040. Each Hne is spHt due 
to the Kepler rotation of the disk but shows an asymmetry in the relative strengths of the red/violet peaks, 
favoring the red wing peak in this case. 



with time, we later modified our spectrum- simulation code according to these geometri- 
cal situations. 



For the hydrodynamic investigation of a viscous gas-disk we used the code FARGO 
(Fast Advection in 7?otating Gaseous Objects) by Masset [9]. Even though being almost 
exclusively used in the calculation of protoplanetary disks and planet formation scenar- 
ios, this code is suited for all kinds of sheared fluid disks. We calculated the temporal 
evolution of a blob of material, starting with a Gaussian density distribution (parameters 
given in Tab. 1). 

To perform the modifications on the spectrum synthesis code, we selected two ge- 
ometries, a spiral arm structure of the early simulation phase and a fully closed but 



Hydrodynamic Simulation with FARGO 



TABLE 1. FARGO simulation parameters 



White dwarf mass Mwd/M© 0.77 

White dwarf radius ^wd /km 7700 

Disk mass Mdisk/g 7-10^' 

Disk outer radius Rq/Rwd 136 

Simulation rim Rsim/Ro 1-5 
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FIGURE 3. The two selected evolution phases of the FARGO simulation - color-coded is the surface 
mass density on a logarithmic scale. The top panel shows an early stage (fsim ~ 40 yr) as the disks forms a 
spiral-arm like structure. The lower panel shows a late {tsim ~ 125 yr), but with several hundreds of orbits, 
long-lasting phase of a fully closed disk with an out-of-balance matter distribution. 



out-of-balance disk from the late stage of the evolution (Fig. 3). 



Spectral Synthesis Code AcDc 



Using the Accretion Disk code (AcDc) developed by Nagel et al. [11] to calculate 
the disks' spectra, we assumed a thin a-disk as described by Shakura and Sunyaev [13]. 
This allows us to decouple the vertical and radial structure and separate the disk into 
concentric annuli. For each of these rings with radius R the effective temperature T^ff as 
well as a viscosity times surface mass density wL value is derived for given values of 
the white dwarf's radius RwD, mass Mwd, accretion rate M, the gravitational constant 
G, and the Stefan-Boltzmann constant a from 



The program then simultaneously solves the following set of equations for each of the 
rings: 

• Radiation transfer for the specific intensity / (v, ju), 

• Hydrostatic equilibrium between gravitation, gas and radiation pressure, 

• Energy conservation for viscously generated £'mech and radiative Erad and 

• Static rate equations ^ = for the population numbers of the atomic level / of 
the given model atom in the non-LTE case. 

For the subsequent spectral surface integration phase of AcDc, we constructed simple 
geometry maps (Fig. 4) following the selected FARGO pictures. This modification al- 
lowed us to reset the spectral flux to zero for those ring segments, which are not part 
of the structure. By simple linear interpolation between the neighboring annuli spectra, 
AcDc integrates the spectrum of the whole non axis-symmetrical disk. 



Calculated with the model atoms given in Tab. 2 and a set of parameters (Tab. 3), the 
spectra of the two geometry maps are shown in Fig. 5. 

As a first result of our analysis, we had to increase the inner radius of the disk to at 
least i?i>58i?wD in order to suppress broad line wings and to fit the flat continuum 
between the triplet components as well as the line's steep profiles. The latest result 

TABLE 2. Model atoms and element abundances 

Ions Hi- II Ci-iv Oii-iv Mgi-iii Sii-iv Cai-iV 

NLTE levels 11 45 32 57 40 38 

Line transitions 45 88 65 104 49 64 

Abundance (%, in mass frac.) 10"^ 4.6 65.5 13.5 15.1 1.3 




(1) 



(2) 



SYNTHETIC SPECTRA AND RESULTS 




FIGURE 4. Geometxy maps representing the selected FARGO phases shown in Fig. 3. The two colors 
(greyshades) indicate that the flux of the adjoint ring segment is set to zero and AcDc will compute a 
downward interpolation. For the subsequent spectrum calculations the observer is meant to look from the 
lower end of the box towards the central object. 
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FIGURE 5. Synthetic spectra for both geometries (shown in Fig. 4), in comparison with the observation. 
The line asymmetry as well as the flat continuum between the lines are fitted well for an inclination of 
i = 76°. The marked C II feature seen to the far right of the modeled spectrum is actually a yet unspUt 
doublet line. 



TABLE 3. AcDc integration parameters 



Ring radii range R/Ryju 
Disk surface mass density (g/ cm-^ ) 
Disk temperature range Tiff/K 



2 - 136 
0.3 

6700 - 5600 



by Melis et al. [10] being R[ — 40±3i?wD is only slightly smaller than our value. As 
matching the asymmetric line profile was the main goal of this work, we were able to 
achieve this well for both geometries, although the relative strengths of the red/violet 
parts of each line seem to be generally estimated too low in the case of the closed disk 
structure. 

Of course the effect of the geometrical asymmetry strongly depends on the orientation 
towards the observer. With the rotation of the asymmetrical disk around the central 
object, the line profile should vary. Gansicke et al. [1, priv. comm.] did, in fact, report 
on a change in the line profile in the order of years for at least two of the known five 
objects. But the predicted variability timescale of our model seems to be much shorter 
(order of hours) than the observations do suggest. 

On the other hand, there might also be a long-time variability as the FARGO simula- 
tions suggest a disk evolution through different geometries in the time of hundreds of 
orbits. 
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